Thermodynamic properties of the one-dimensional Kondo insulators 
studied by the density matrix renormalization group method 
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Thermodynamic properties of the one-dimensional Kondo lattice model at half-filling are studied 
by the density matrix renormalization group method applied to the quantum transfer matrix. Spin 
susceptibility, charge susceptibility, and specific heat are calculated down to T = O.lf for various 
exchange constants. The obtained results clearly show crossover behavior from the high temperature 
regime of nearly independent localized spins and conduction electrons to the low temperature regime 
where the two degrees of freedom couple strongly. The low temperature energy scales of the charge 
and spin susceptibilities are determined and shown to be equal to the quasiparticle gap and the spin 
gap, respectively, for weak exchange couplings. 
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The Kondo lattice model (KLM) is a simple theoretical 
model for heavy Fermions which consists of two different 
types of electrons; the localized spins whose charge de- 
grees of freedom are completely suppressed, and the con- 
duction electrons that propagate as carriers in extended 
orbitals. Both the conduction electrons and the local- 
ized spins do not interact among themselves, however, 
the exchange interaction between them leads to correla- 
tions and yields various interesting physical phenomena. 

Recent studies on the one-dimensional KLM has shown 
that this model is always insulating at half-filling, [y This 
insulating phase has different excitation gaps for the spin 
and charge channels. The spin gap defines the energy cost 
needed for the lowest spin excitation which changes the 
total spin quantum number. The charge gap defines that 
for the lowest pure charge excitation which changes the 
total carrier number by two keeping the spin quantum 
numbers fixed. 

The spin gap of this insulating phase has been exten- 
sively studied both analytically and numerically. ^^ 
It has been shown that the spin gap Ag has similar J- 
dependence to that of the Kondo temperature Tk for 
the single impurity Kondo model. The numerical studies 
show As oc exp(— 1/apJ), where p is the density of state 
at the fermi level and a is the lattice enhancement fac- 
tor. The lattice enhancement factor was recently deter- 
mined to be a = 1.4 by the density matrix renormaliza- 
tion group (DMRG) method. Q The similarity between 
the Kondo temperature and the spin gap suggests that 
the origin of the spin gap is due to some sort of singlet for- 
mation between the conduction electrons and the local- 
ized spins which is the essence of the Kondo effect. The 
difference appears in the coefficient a in the exponent, 
which is due to collective singlet formation of conduction 
electrons with many localized spins in the lattice. 

For the charge gap, on the other hand, linear J- 
dependence (Ac = J/2) has been obtained for weak ex- 
change couplings. MM This J-dependence may be under- 



stood by the strong antiferromagnetic correlations among 
localized spins that generate a staggering internal mag- 
netic field on conduction electrons which is proportional 
to the exchange constant. It should be noted that the 
correlation length of the spin degrees of freedom is much 
longer than the charge correlation length. Therefore the 
above arguments are justified in spite of the fact that 
there is no magnetic long range order. 

The qualitatively different nature of the spin and 
charge gaps is a consequence of quantum mechanical 
many body effects among the conduction electrons and 
localized spins, and is a unique feature of this strongly 
correlated insulating phase. 

In the present paper we will study how such inter- 
play between the conduction electrons and the localized 
spins appears in the temperature dependence of various 
thermodynamic quantities. To discuss the thermody- 
namic properties we use the finite temperature density 
matrix renormalization group method (finite-T DMRG). 
We calculate spin and charge susceptibilities and spe- 
cific heat for various exchange constants. We will show 
clear crossover behavior between high and low tempera- 
ture regimes. In the high temperature regimes the local- 
ized spins and the conduction electrons are only weakly 
coupled, while in the low temperature regime they couple 
strongly. Characteristic energy scales in the low temper- 
ature regimes are identified as the spin gap and the quasi- 
particle gap for the spin and charge degrees of freedom, 
respectively. 

The Hamiltonian we use here is the one-dimensional 
KLM described by 



H 
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where the operator Cis (cj^) annihilates (creates) a con- 
duction electron at site i with spin s (=T,i) {S^^ = 

(1/2) X^s s' '^Is^ss''^is')- The hopping matrix element is 
given by t and J is the antiferromagnetic exchange cou- 



pling between the conduction electron spins Sc and local- 
ized spins Sf, both beeing spin 1/2 degrees of freedom. 
The density of conduction electrons is unity at half- filling. 

In order to calculate thermodynamic quantities we use 
the density matrix renormalization group method [pj ap- 
plied to the quantum transfer matrix. Recently this 
method of finite-T DMRG was successfully applied to 
the one-dimensional quantum spin systems to calculate 
thermodynamic quantities. |^,|l^ The present study is 
the first application of the finite-T DMRG to a system 
with Fermionic degree of freedom. 

In the following we briefly outline the method. We 
define the transfer matrix 
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where M is the Trotter number. Here the Hamiltonian H 
is decomposed into two parts i?odd = X]n=i f^2n-i,2n and 

Heven = z2n=l ^"^^,271+1 SUCh aS [/l2n-l,2n, ^2n'-l,2n'] — 

[h2n,2n+i, /i2ra',2n'+i] = 0. First we diagonalize the trans- 
fer matrix with a small M to obtain the maximum eigen- 
value A and corresponding eigenvector. Next we calculate 
generalized density matrix from the obtained eigenvector. 
By diagonalizing the density matrix, we chose important 
basis states which have large eigenvalues for the repre- 
sentation of the transfer matrix. Using these basis states 
we increase the Trotter number M of the transfer matrix 
within the fixed number of basis states and continue the 
above procedure until we get the Trotter number suffi- 
cient for a given temperature T. The free energy of the 
infinite system is directly obtained from the maximum 
eigenvalue A of the transfer matrix: F = — (T/2)lnA. 
The spin and charge susceptibilities, and specific heat 
are obtained by numerical derivatives of the free energy. 

The following calculations are performed by the infi- 
nite system algorithm of the finite-T DMRG keeping 40 
states per blocks. [g|,|0[ The truncation error in the finite- 
T DMRG calculation is typically 10"^ and 10"^ at the 
lowest temperature with the Trotter number M — 50. To 
check the results we have compared the obtained spin and 
charge susceptibilities to those calculated by the Quan- 
tum Monte Carlo simulations for J/t = 1.6. |0] The 
overall structure agrees well, but low temperature part of 
our results are more reliable. For example, the estimated 
spin gap energy 0.7t in Ref. 11 is not consistent with 0.4i 
obtained by the standard zero-temperature DMRG [g, 
while it is consistent with the present result (0.45 ± 0.1)i 
obtained by the finite-T DMRG. 

We first show the temperature dependence of the uni- 
form spin susceptibility. The spin susceptibility is ob- 
tained from the change of the free energy by a small 
magnetic field h: SF — Xs/j^/2. The results for J/t — 
0, 1.0, 1.2, 1.6 and 2.4 are shown in Fig. 1. In the limit 
of J/t = 0, both the localized spins and conduction elec- 
trons are uncorrelated. The susceptibility is given by the 
sum of the Curie term of localized spins and the Pauli 
term of free conduction electrons. In this system the con- 
tribution of the Pauli susceptibility of the free electrons 



(shown in Fig. 1) is relatively small, and the total sus- 
ceptibility of J/t = is dominated by the Curie term of 
the localized spins, which diverges in the limit of T = 0. 
With introducing the exchange coupling, low temper- 
ature part of Xs sharply decreases with decreasing the 
temperature. This drastic change is due to the forma- 
tion of the spin singlet state between the localized spins 
and conduction electrons whose energy scale is given by 
the spin gap for small exchange constant. The spin gap 
As obtained by the zero temperature DMRG is 0.08t for 
J = l.Oi, which is consistent with the characteristic tem- 
perature at which Xs starts to decrease. 



?>J 



1.4 


i' ' ' ' ' ' 


- 


1.2 


\ 


- 


1 


j/t=\.a\ \ 


- 


0.8 


r\\\ 




0.6 


/0/r=li'^N\. \^^ 


- 


0.4 
07 


.' "-^^^^^tl^-,.,,,^^^ "^■~'- „.^ free fermion 


^ 











: .-I--"' 1 1 1 1 1 


1 







0.2 



0.4 



1 



1.2 



1.4 



0.6 0.8 

T/t 
FIG. 1. Spin susceptibility of the half-filled 
one- dimensional Kondo lattice model. The truncation error 
in the finite-T DMRG calculation is typically 10"^ and 10"^ 
at the lowest temperature. 

This characteristic temperature separates high temper- 
ature region where Xs increases with decreasing the tem- 
perature, and low temperature region where the suscep- 
tibility drops rapidly with decreasing temperature. This 
temperature is expected to be determined by the com- 
petition between the thermal fiuctuations and the sin- 
glet correlations. Since the exchange interaction stabi- 
lize the singlet correlation between the localized spins 
and conduction electrons, the crossover temperature in- 
creases with enhanced exchange coupling. 

In order to determine the energy scale at low tempera- 
tures more precisely we estimate the activation energy by 
fitting the obtained susceptibility. The estimated activa- 
tion energy for the spin susceptibility is summarized in 
Table I for J/t — l.Qt,2At, and 3.0i. Here we can com- 
pare this energy scale with quasiparticle gap and spin 
gap, both responsible for the magnetic excitation. The 
quasiparticle gap and spin gap are obtained by the stan- 
dard zero temperature DMRG method. 

In Table I we see that the lower one of the quasiparticle 
gap and the spin gap determines the low temperature 
energy scale of the spin susceptibility. This is consistent 
with the general form of the susceptibility that is formally 
written as 
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According to the above equation the low temperature 
behavior is determined by the energy difference between 
the ground state and the lowest excited state which is 
active for magnetic field. 

Next we consider the charge susceptibility. The charge 
susceptibility Xc is obtained from the change of free en- 
ergy by a small shift of chemical potential /i, 5F = 
Xc\^ 1^- In the present calculation we use the fact that 
the chemical potential is zero at half-filling owing to the 
SO (4) symmetry of the model. 

Calculated Xc for Jjt = 0,1.0,1.2,1.6 and 2.4 are 
shown in Fig. 2. For J/t = 0, Xc does not show diverging 
behavior at low temperatures in contrast to Xs- In the 
limit of T = Xc is equal to the density of state of con- 
duction electrons that is I/tt. This is natural since the 
charge degrees of freedom is solely governed by the con- 
duction electrons. Since there is no correlation between 
up-spin and down-spin conduction electrons for J/t = 0, 
Xc/4 is equal to the spin susceptibility of the free con- 
duction electrons. The slight increase in Xc in the low 
temperature region is a characteristic feature of the one- 
dimensional system where the density of states diverges 
at the band edge. 
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FIG. 2. Charge susceptibility of the half-filled 
one-dimensional Kondo lattice model. 

Switching on the exchange coupling, Xc decreases 
rapidly at low temperatures. Similarly to Xs, this is due 
to the singlet formation between the conduction electrons 
and the localized spins. With increasing the exchange 
coupling, decrease in Xc appears at higher temperature. 
This is consistent with the behavior of the charge gap 
(quasiparticle gap), which is enhanced with increasing 
the exchange coupling. 

To determine the energy scale at low temperatures. 



we estimate activation energy for the charge suscepti- 
bility. By fitting Xc with an exponential function, ac- 
tivation energy A^^ is obtained as listed in Table I for 
J/t — 1.6t, 2.4t, and 3.0t. From this table it is concluded 
that the quasiparticle gap always determines the low tem- 
perature energy scale of the charge susceptibility. Since 
the charge gap is twice the quasiparticle gap, the lowest 
excitation that is responsible for the charge excitations 
is always the quasiparticle gap. 

Finally, we consider the specific heat. The specific heat 
is calculated from the second derivative of the free energy; 
C = -Td^F/dT'^. The resuhs for J/t = 0, 1.0, 1.2, 1.6 
and 2.4 are shown in Fig. 3. 

At J/t — the specific heat of this model is given by 
the sum of the delta function at T = that originates 
from the localized spins and the specific heat of free con- 
duction electrons. By including the exchange coupling, 
they are combined to make a two-peak structure. The 
peak at higher temperatures is almost independent of 
the exchange constant, and similar to the structure of the 
free conduction electrons. Thus, this structure originates 
from the band structure effect of the conduction elec- 
trons. The peak at lower temperatures strongly depends 
on the exchange constant. The peak shifts toward higher 
temperatures and becomes broader with increasing J/t. 
This behavior is consistent with that of the spin suscepti- 
bility whose peak also shifts toward higher temperatures 
and becomes broader with increasing J/t. Thus, we can 
assign the origin of the peak at lower temperatures due 
to the localized spin degrees of freedom coupled through 
the conduction electron spins. The slight decrease in the 
specific heat at low temperatures for weak J/t compared 
with that of the free conduction electrons is due to the 
opening of the charge gap. The decrease at low tempera- 
tures are compensated by the increase at higher temper- 
atures as shown in Fig. 3. 

Further increasing the exchange coupling, large exci- 
tation gap which is comparable to the hopping matrix 
element t opens for both spin and charge sectors. The 
ground state is close to the collection of the local sin- 
glets, and the specific heat has a single peak structure. 
(See J/t=2.4 in Fig. 3) 

TABLE I. Activation energy obtained from the spin and 
charge susceptibility, A^^ and A;^^ , and the quasiparticle gap 
Aqp and the spin gap As of the one-dimensional Kondo lat- 
tice model. The charge gap is twice the quasiparticle gap; 
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FIG. 3. Specific heat of the half-filled one-dimensional 
Kondo lattice model. 

In conclusion we have successfully applied the finite- 
T DMRG to the one-dimensional Kondo lattice model 
at half-filling. Temperature dependence of the spin and 
charge susceptibility as well as the specific heat are cal- 
culated down to T = O.li. Compared with the quantum 
Monte Carlo simulations, an advantage of the finitc-T 
DMRG is that it is free from statistical errors and nu- 
merical errors in the calculations are under control to cer- 
tain extent by keeping sufficient number of basis states. 
More importantly it is, in principle, free from the nega- 
tive sign problem which practically invalidates the quan- 
tum Monte Carlo simulations for frustrated spin systems 
and most of Fermionic systems. 

The subtle interplay between the localized spins and 
the conduction electrons are clearly seen in the temper- 
ature dependence of the susceptibilities and the specific 
heat. The low-temperature energy scale of the spin sus- 
ceptibility is determined by the smaller one of the spin 
gap and the quasiparticle gap. Thus, the spin gap is 
the low energy scale in the weak coupling regime. On 
the other hand, the low energy scale for the charge sus- 
ceptibility is always determined by the quasiparticle gap, 
which is half of the charge gap. Effects of both the spin 
gap and the charge gap are seen in the specific heat. 

The present study is the first one where the finite-T 
DMRG is applied to a system with Fermions. Existence 
of the excitation gap at half-filling is favorable for con- 
vergence of numerical calculations. Of course it is an 
interesting future problem to apply the finite-T DMRG 
to the Kondo lattice model away half-filling where the 
ground state is a Tomonaga- Lut t inger liquid with gap- 
less spin and charge excitations |l^ , |l3| . 
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